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CrossMark 
Abstract 
We investigate the wave propagation on a compact 3-manifold of constant 
positive curvature with a non trivial topology, the Poincaré dodecahedral 
space, if the scale factor is exponentially increasing. We prove the existence 
of a limit state as t > +00 and we get its analytic expression. The deep sky is 
described by this asymptotic profile thanks to the Sachs—Wolfe formula. We 
transform the Cauchy problem into a mixed problem posed on a fundamental 
domain determined by the quaternionic calculus. We perform an accurate 
scheme of computation: we employ a variational method using a space of 
second order finite elements that is invariant under the action of the binary 
icosahedral group. 
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(Some figures may appear in colour only in the online journal) 


1. Introduction 


The topology and the geometry of our universe is a fundamental open problem. There are a lot 
of articles about these topics (for example [23, 24, 35]). Recently the data from the WMAP 
(Wilkinson Microwave Anisotropy Probe) mission have provided strong evidence suggesting 
that the observable universe is nearly flat, with the ratio of its total matter-energy density to the 
critical value very close to one [7], but without fixing the sign of its curvature. So the study of 
the wave propagation in universes involving the Poincaré dodecahedral space continues to be 
relevant in particular for the search of the circles-in-the-sky signature [1, 2, 27]. This model 
has received a lot of interest and generated numerous works (see e.g. [1, 5, 8, 20, 25, 30, 36]). 
These spacetimes are Lorentzian manifolds (R, x K, g) that are globally hyperbolic and multi- 
connected. Here K is the Poincaré dodecahedral space that is a three-dimensional C% compact 
manifold, without boundary, defined as the quotient of the 3-sphere S? under the action of the 
binary icosahedral group Z* ([14, 31, 34, 37]). In [5] we have investigated the propagation of 
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the scalar waves in the case of the stationary universe for which ds? = dt? — dsk where dsk 
is the Euclidean metric induced by the three-sphere S? on K. In this paper we study the wave 
propagation on dynamical universes (U, g) of Robertson—Walker type 


U=R,xK, dsy = g,,dx!'dx” = dt? — a?(t)dsk (1) 


where the scale factor a(t) is a smooth positive function. We are mainly interested in consider- 
ing the accelerating spacetimes for which the scale factor is an increasing function with an 
exponential behaviour, a(t) ~ e% as t > oo. In the study of the asymptotics of the fields at the 
infinity, and for the numerical experiments, we treat two important cases. In the first one, we 
choose a(t) = H~'cosh(Ht) where H > 0 is the Hubble constant. Then g is solution of the 
vacuum Einstein equations with cosmological constant A = 3H? and U is locally isometric 
to the de Sitter spacetime dS, and differs globally from dS, by its topology that is not trivial. 
We denote this universe dS4(K). In particular, while dS4 is spatially homogeneous and glob- 
ally isotropic, dS4(K) is homogeneous but not globally isotropic [26]. In the second model 
we choose a(t) = e*. Then g is no longer a solution of the vacuum Einstein equations but it 
is a toy model for a universe with an inflation generated by a scalar field or a solution in f(T) 
theory (see e.g. [16]). 

In this paper we investigate the scalar waves on (U/,g) that are the solutions of the 
D’ Alembert equation 


1 
OY := a(o Jel av) =0. 
lel 


For the spacetime (1), this equation has the form 


at) 5 1 
= 


E +3 == 
a(t) a(t) 


a| Y= 0, (2) 
where Ax is the Laplace-Beltrami operator on K. This equation plays a fundamental role in 
the theory of cosmological perturbations (see e.g. [28]). It arises in the linearized Einstein 
equations in the vacuum, written in harmonic coordinates, and in this case W is the fluctua- 
tion around the ground metric, i.e. Y describes the gravitational waves. On the other hand the 
Sachs—Wolfe theorem assures that in the Friedmann—Robertson—Walker (F—R—W ) universe 
with an ultrarelativistic gas, the scalar perturbations obey to the acoustic equation. This result 
was proved originally for the spatially flat metric, and has been extended to arbitrary space 
curvature [12, 15]. In this context, considered in [1, 2], we write the ground metric (1) as 


ds? = A(n) [dn? — dsk], dsk = dy? + sin? y(d6? + sin? Ody’), 


where n = fi a` }(t)dt is the conformal time, and A(n) = a(t) is the cosmic scale factor. Then, 
the metric with the scalar perturbation Y takes the form in conformal Newtonian gauge: 


ds? = A(n) [C1 + 2W)dr? — (1 — 2W)dsx) . (3) 
The knowledge of Y allows one to adress a fundamental issue of the cosmology: the study 
of the temperature fluctuations 67 (6, p) of the temperature T(0, y) of the cosmic microwave 
background radiation, observed in the direction (0, Y) € S°. The integrated Sachs—Wolfe form- 


ula (see i.e. [1]) gives a link between 67/T observed at the conformal time ng of our present 
epoch, and Ų between ng and n the conformal time at recombination ([1], formula (31)). 
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Therefore we must know the time evolution of Y to be able to use it. An approximated formula 
for low spherical harmonics is also given by the ordinary Sachs—Wolfe formula, 
éT(6, 9) I 


~ 3 Pe X = No — Mis 8, P). (4) 


If m- > arecos( "= ) where ø is the golden ratio, the horizon radius of the observer is 
larger than the injectivity radius, hence repeating images of the sphere horizon intersect each 
other. Since the self-intersections of the sphere of last scattering are six pairs of matched antip- 
odal circles, the fluctuations of a realistic Y(n; X = No — Ns -), for instance if the primordial 
fluctuation amplitude (0, .) is a Gaussian random variable, would be correlated around these 
circles: the observer can see circles-in-the-sky [11]. 

Besides, we remark that, as regards the propagation of the electromagnetic fields in the 
curved space-times, the D’Alembert operator O = V°O, is also the principal part of the 
Maxwell equation for the vector potential V°V,A? — RSA = J°. In concluding, there are 
numerous motivations to study this equation, represent and compute the solutions, investigate 
their properties. 

The explicit form of the solutions is known for the classic F-R—W space-time with the 
R x S° topology [19]. The case of the dodecahedral space is much more complicated. In this 
work we solve the Cauchy problem, we represent the solutions in terms of expansions involv- 
ing the eigenmodes of the Laplacian on K, then we investigate the asymptotic behaviour of 
W at the time infinity, and finally we develop a numerical scheme of accurate computation. A 
crucial result is the exponentially fast convergence of V(t, X) to a final asymptotic state V(X) 
when ż tends to the infinity. The link between Y and the initial data of Y is given by a pseudo- 
differential operator. In particular for dS4(K), we have 


T= ee bs sin( = JD F TH, J+ a Í. ƏT, Jwg, (5) 


where wxg is the volume form on K. To compute Y with this formula we should use the eigen- 
functions of the Laplacian. Despite significant progress for the determination of these modes 
(see in particular [6, 20]), the implementation of this strategy seems to be challenging. Instead, 
we solve numerically the initial value problem, and taking advantage of the fast convergence, 
we simply find the asymptotic profile as 


U(X) = lim Wt,X)~+W(T,X), T> 1. 
t> +00 


From the point of view of the theory of cosmological perturbations, this result has a great inter- 
est if Y is the scalar perturbation in the metric (3). If the conformal time of recombination is 
large enough, we have Y(n, .) ~ Y. Then the approximated Sachs—Wolfe formula becomes 


6T(@, 9) _ I 
T 3 
and the circles-in-the-sky have simply to be found in Yy- 
We briefly describe the following sections. In part 2 we solve the global Cauchy problem 
for (2) in the functional framework associated to the energy 
1 
a(t) 


W(x = No a Nis» 0, p), (6) 


EY, t) = 3I Olit DAORA (7) 
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Figure 1. The Fundamental Domain (dark blue). We display F, (in dark blue) and the 
projection p in R? of the images of F by six elements of Z*. It is the representation of a 
part, included in B’, of the tiling of S°. 


where the covariant nabla operator satisfies < VkY, Vk® >p) =< —AkY, ® > pw) 
for any Y, € C™(K). We represent the Poincaré dodecahedral space by its fundamen- 
tal domain F that is the geodesic convex hull of twenty vertices on S° that are included 


2 
in {xo = ap x R? where ø is the golden ratio. We introduce the domain of vizualisation 


F, C R? that is the projection (xo, x1, X2, 23) E S 31+ (xp, X3) € R? of the fundamental domain. 
F, is a dodecahedron with pentagonal curved faces included in ellipsoids (figure 1). Defining 
u(t, Xi, %2, X3) := WF, xo, X1, %2, X3) the wave equation (2) is equivalent to 


‘ 1 
a+ 3400 Ag|u=0, tER, xEF, 
[o a) ao y A 


together with the boundaries conditions 
YTET*, (X0 Xp X2 X3) = T(X0, X, X, X3) => U(t, X1, X2, X3) = U(t, X1, 32X3). (9) 


We solve the initial value problem associated to (8) and (9) and we state an equivalent varia- 
tional problem. We discuss the existence of a future horizon and we compute the radius of the 
horizon sphere. 

In the third section we investigate the asymptotic behaviour of W(t, .) as t > +00. We prove 
the existence of an asymptotic profile Y% := lim,_,,, Y(t) if the scale factor is exponentially 
increasing, and we compute an analytic expression of V,, in the cases a(t) = H~! cosh(Hr) and 
a(t) = e. These results extend those of [3] regarding the steady state universe (see also [33]). 
We estimate the norms of W,, by the energy of Y. 

In the fourth part we numerically investigate the mixed problem (8, 9). We are mainly 
interested in computing the waves for a long time, hence the exponential behaviour of the 
coefficients in the equation leads to a hard numerical challenge. To overcome this difficulty we 
have to perform very precise computations. To reach this goal we adopt the variational method 
associated to (8) with suitable finite elements satisfying constraint (9). In the case of the sta- 
tionary universes we used in [5] finite elements of P, type. To be able to treat the accelerating 
universes, we need more precision and we use elements of second order. Our numerical results 
are in excellent agreement with the theoretical results on the radius of the future horizon and 
on the asymptotic behaviours of the waves. 
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2. Waves on dynamical universes 


In this section we solve the global Cauchy problem for the D’Alembert equation on the 
dynamical universes (1). We assume the scale factor a is aC™ function that is strictly positive. 
As a direct consequence of the theorem of Choquet-Bruhat, Cotsakis (theorem 11.10 of [9]), 
we can see that (U, g) is a globally hyperbolic manifold. Therefore according to Leray [22], 
the global Cauchy problem is well posed in C%. Now it is natural to look for the solutions in 
the scale of the Hilbert spaces associated to the energy (7). Given m € N, we introduce the 
Sobolev space 


H™K) := {Y € L(K), VẸY € LK), Jal <m} 


endowed with the norm 


ltir: = YS IVE Za. 


ja|<m 


Here Vx are the covariant derivatives and L?(K) is the usual L?-space of Lebesgue associated 
to the volume form on K. We can also interpret this space as the set of the distributions W that 
belong to the usual Sobolev space H” (S°) such that Y o 7 = Y for any 7 € Z*. We use also the 
space H~!(K) that is the dual space of H'(K). As usual, we identify Z? with its dual space and 
we have H'(K) C L(K) C H~'(K). The subspace C™(K) is dense in all these spaces and we 
have for $, Y € C% 


(Y, ®) y- 1K), HK) — = (4*, Pew 
We also know that —Ax + 1is an isometry from H! (K) onto H- (K). We shall use an orthonor- 
mal basis in Z?(K), formed of eigenfunctions W} C H! (K) of the densely defined self-adjoint 
operator — Ax associated to the eigenvalues q’, i.e. -AKW, = q?W,. If W(X) = gAqld(X), 
we have 

Illia = DQ? + D"A P. 

q 
A {finite energy solution’ is a solution Y of (2) that belongs to C°(R,; H'(K)) N CIR; L’(K)). 

Since Ax is continuous from H! to H~!, we remark that sucha solution belongs to CR; H (K)) 


since av = = pbk 320 aW) € CHR; HK). 


Theorem 2.1. Given m €N, Yo € H” (K), Yı € H”(K), t* € R, there exists a unique solu- 
tion Y of the wave equation (2) satisfying 


Y € C'R; A"*(K)) N C'R; AK), (10) 
Vr, = W), OM.) = WA. (11) 


For m = 0, Y is the unique solution satisfying together (10), (11), and for all ® € H'\(K) 


a(t) d 


W(t), ®) 2 
x (t), ®) pay + 3—— at) di 


TOO. Dr + 2 Tp KUO. Veð) 2a = 0. (12) 


If the scale factor a(t) is increasing on (t*, o0), then the energy defined by (7) is decreasing 
on this interval, and if a(t) also tends to the infinity as t > œ, then the energy tends to zero. 
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Remarks. If m=0 we have C?(R,;H~'(K)), hence for any ®€A'(K) we have 
(W(t), ®) pay = (Y(t), ®) 149, ma E CR). Therefore the first term in (12) is well 
defined. If Yo, Yı € C™(K), then YE C%™(R, x K) by the Sobolev embedding theorem 


m H”(K) = C™(K), and we say that Y is a ‘smooth solution’. If the energy tends to zero, 
we know that ||0;Y(4)||72qx) tends to zero, but we have no control on || Y(4)||;2¢q). We shall prove 
in the next section that W(r) tends to an asymptotic state UV, and the norms ||W,ol|;2¢q) and 
|| Vi Yollz2aq) can be estimated by the norms of the initial data. 


Proof. There are a lot of ways to obtain the theorem from well known results. To construct 
the solutions we could invoke the existence of the smooth solutions by the theorem of Leray 
[22], then we extend by density of C% in H’(K) by using the energy estimate satisfied by the 
regular solutions: 


a'(t) 


a(t) 


d t 
Tkw = 6f n NOVO) Begg — 2 IV WOO) [Bagg (13) 


a 


Then by the Grönwall lemma we can control the energy at any time, hence we get the global solu- 
tions for m = 0 (see also theorem 2.29 in appendix III of [9]). Finally commuting the covariant 
derivatives and the D’Alembert operator, we deduce the existence of solutions for any m € N. 
Another way consists in employing the spectral method of Kato (see e.g. [32]). In fact the more 
convenient approach to directly obtain the existence, the uniqueness and the equivalence with the 
variational problem is the variational method by Lions.There are obtained by a direct application 
of the results of chapter XVIII of [13]. Finally if a’ is positive, (13) implies that for all t > t* 


o E(W,t), 


I kij- 
dt alt 


hence the energy is decreasing and “(a(EW, t)) < 0. We deduce that 


E(Y,t)< a) row, t*). 
a(t) 


We conclude that the energy tends to zero if a(t) > oo. oO 


To perform a computational method to solve the Cauchy problem, it is convenient to rep- 
resent the Poincaré dodecahedron K by a solid £F, in R?. We recall that K is the quotient of the 
3-sphere S? under the action of the discrete fixed-point free subgroup Z* of isometries of S?, 
with Z* acting by left multiplication (see [1, 5, 8, 20, 25, 30]). The subgroup Z* of SO(4) is 
the binary icosahedral group ([14, 31, 34, 37]). It is a two sheeted covering of the icosahedral 
group Z C SO(3) consisting of all orientation-preserving symmetries of a regular icosahedron. 
To visualize S? we use the parametrisation: 

Xo = COS X 

xı = sin y sind sin Y 

X% = sin x sin 0 cos yp 

x3 = sin x cos 6. (14) 


Hence S? can be seen as two unit balls, B’, B” in R? glued together by their boundary: 


S? = (B'U B"IS”. (15) 


Class. Quantum Grav. 34 (2017) 055010 A Bachelot-Motet and A Bachelot 


The radial coordinate in these balls is given by r = sin x. The 0 and ọ coordinates are the 
standard ones for S*. For the first ball x runs from 0 at the center through 5 at the surface. For 
the second ball y runs from = at the surface through 7 at the center. 

To construct F,, we represent S*/Z* by a fundamental domain F C S? and an equivalence 
relation ~ such that 


Sl = Fix, (16) 
The construction is based on the quaternionic calculus. It is detailled in [5]. Here we just pres- 
ent the main features. F is such that: 


s=U+@), and YvreT* Vr'e€l*, rer sl(PDNMD =Ø. (17) 
TET* 


F is a regular spherical dodecahedron (dual of a regular icosahedron). One hundred twenty 
such spherical dodecahedra tile the 3-sphere in the pattern of a regular 120-cell. More spe- 


cifically we choose the fundamental domain such that (1,0,0,0) € F, then F is the geo- 


desic convex hull in S? of these 20 vertices: (0, -,0 1), (0, 0,+1, 2y 


o? 


2 7 (a°, = F A and a 
Each of the 12 faces of F is a regular pentagon in S? (see [5], for details). 

The equivalence relation ~ is defined by specifying the equivalence class q of q € F. It is 
built from twelve Cliffort translations g; that have been used for the construction of F (see the 


appendix and [5]): 
å = (g(9), i€ {1,..., 12}}F. 


It follows that: (1) if g belongs to F , then ġ has only one element, (2) if q is a vertex of F, 
then ġ has four elements, (3) if q belongs to an edge of a face, without being a vertex, then g 
has three elements, (4) if q belongs to a face and does not belong to an edge, then g has two 
elements. 

Finally to define F, we use the projection 


(o?, +1, = 0), where o = (1 + V5) /2 is the golden number. 


P+ (X0,%1, 20, 13) € SS — (241,29, x3) € RY, (18) 
and we choose 
Fy := pF). (19) 


With our choice of fundamental domain p is one-to-one on F and we have 


x.y, 2) ER (Jl - x- y? — 2? x,y, EF. (20) 


In figure 1, F, is depicted in dark blue, as a part of B’. The vertices of F, are those of a cen- 


tered regular dodecahedron in R? and they are given by =a 5, 0,1), 55 0. 1, >, 


z L, L, *) and xt 1, 5,0). Each face of F, is included in an ellipsoid [5]. The 


equivalence relation ~ on F induces canonically an equivalence relation ~ on Fy: 


VX, X' EF, X= XS p(X) ~ p(X’). (21) 


~ leaves invariant the interior of F, and identifies any pentagonal face of F, with its opposite 


face, after rotating by 5 clockwise around the outgoing normal at the center of this last face. 
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Obviously W is entirely determined on R,xK_ by its restriction to R,x F, 
hence we are going to work on R, x F,. We introduce the map f: (x,y,z) E Fu f(x, y,z) = 


(J l= y —2’, x,y, z) € S?. fis one-to-one from F, C R? onto F C S? and we put 


ult, x,y,z) = V(t, f(x, y, z)). (22) 
It is clear that 6 € H'(K) iff v := ® o f belongs to the usual Sobolev space H! on F, for the 
euclidean metric of R?, H(F) = {v € L(A), Vey Vv € L’(F,)}, and satisfies the boundary 
condition 

VX,X'ECOF,, X~X! => W(X) = X’). 
We deduce that Y E€ CYR; HK) N CIR; L’(K)) is solution of (2) iff u belongs to 
C°(R,; H\(F,)) O CIR; L’(F,)) and satisfies the equation (8) where 
Ag, = (1 — x?) 62, + (1 — y?) 0% + (1 — 22) 033 — 2xy OF — 2xz OF— 2yz O33 — 3x O1— 3y O2 — 3z Os, 

(23) 


and the boundary conditions 

V(t, X,X') ER x OF, x OF, X~ X'=> u(t, X) = u(t, X’). (24) 
To handle the boundary when applying the finite element method, it is very convenient to take 
into account the boundary condition (24) by a suitable choice of the functional space. In the 


sequel we denote X := (x, y, z), |X|?:= x? + y? + z? and dX := dx dy dz. For m = —1, 0, 1, we 
introduce the spaces W”(F,) that are isometric to the spaces H”(K) : 


W"(F,) = {u= bof, BEH™K)}, lullw":= Nl®lly. (25) 
We have 
WYE) := (Fn (1-|XPY 24X), (26) 
WF) = [uE HF), X~ X uX) = uX), [ul2v= E auoe. (27) 
jaļ<1 


Moreover u € C*(R,, W” (F) iff Y € CR; H”(K)) and theorem 2.1 can be expressed in 
terms of u as the following 


Theorem 2.2. Given uo € W\(F,), u € W(F,), t* € R, there exists a unique u solution of the 
equation (8) and satisfying 


u E€ CORF; WF) N CRI; WAD N COR; WA), (28) 
u(t*,.) = uo(.), ult, .) = m(.). (29) 
u is the unique function satisfying (28), (29) and such that for any ¢ € W\(F,), we have: 


a(t) d 
a(t) dt 


Ë got at 
o= 5 f, G-ixPy tun xnocryax +3 EE f -XP 2u X00 ax 


ta f. AXEI? Vrut, X) WOO = X- Veut, VK Ves AK GO) 
a(t) JF, 


where Vyu = (O,u, Oyu, Ozu), dX = dxdydz. 
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We end this part by discussing the possible horizons in the dynamical universe (1). We 
recall that the comoving spatial distance d between any two points x and y on S? is given by: 
d(x,y) = arccos[6,x'y/]. F, is endowed with the distance induced by d. Moreover the largest 
distance dmax between the origin and the boundary of F, is reached at any of its vertices. So 


2 
o2 


2/2 
cal coordinates p = |X|, w = p7!X, by ds? = dt? — a0 ( dp? + pdu?) where w € S?, 
0 < P< Pmax(w). We have 


4 
Rmax = SUP Prag) = sin(dmax) = ,{1- S ~ 0.378. (31) 
wes? 8 


Given t*ER, O<KR< Rmax we investigate the future causal domain O*(t*,R) of 
{t*} x (AN (0< p < R] x S). Solving 


PO... 1 


J-P att) 62) 


we have dmax = arccos( ja 0.388. The metric on R; x F, can be expressed in spheri- 


we get that 


OH, R) = U {1} x (ANOS pP SRO] x SD) 


where 


R(t) = sin arcsin + f as] (33) 


‘1 ds < 1/2. Therefore an interesting phenomenon appears if 


provided that arcsin(R) + J E 


arcsin(R) + 1 | Aji < arcsin(Rmax) = dmax- (34) 
* als) 


In this case, a future horizon exists: 
Ot(t*,R) C [t*,00) x [0, Ra] x S?, (35) 


where the radius of this future horizon is 
; : co | 
Ry, := sin| arcsin(R) + —ds |]. (36) 
# als) 


Returning to the scalar waves, we conclude that if uo and u; are supported in p < R, then for 
any t > t*, u(t) is supported in p < Rp. 

Similarly, if tı is the time of recombination, there exists a past horizon, called horizon 
sphere or last scattering surface, for an observer located at fobs > tıs. In the universal covering 
S? of K described in spherical coodinates (£, 0, p), itis the submanifold {X = Mobs — Mi} X Sh 3 
where 77 is the conformal time, 


fobs 1 
obs ~ Nis = —— ds < dmax- 
Mobs — Mh J, ie (37) 


Class. Quantum Grav. 34 (2017) 055010 A Bachelot-Motet and A Bachelot 


Then, the horizon sphere on which this observer can see the primordial plasma in the sky S, ie 


i ftos | 
Riptos = SINI f —ds |. (38) 
ts a(s) 


If obs ~ Ms < max we have 


has the radius 


Ritas < Rmax = 4] 1 — — (39) 


fobs 
hence the observable universe of radius R,,, ;,,, 18 strictly included in the whole universe K. In 
this case, there is no multiple copies of images. In contrast, if 


Ri, > Rmax =,/1-— — > (40) 


fobs 
the horizon radius is larger that the injectivity radius. Then the horizon sphere wraps all the 
way around the universe and intersects itself. The case (40) has a great interest in cosmol- 
ogy: the observer could detect multiple images of radiating sources, some circles-in-the-sky 
appear, leading evidence of a non-trivial topology of the Universe [2, 27, 30, 36]. We present 
a numerical experiment of this situation in part 4.2. 

Finally we give the expression of these quantities for the two accelerating universes that we 
consider, dS4(K) and the inflationary spacetime. For a(t) = H~! cosh(Ht), we have 


R(t) = sin(arcsin R — 2 arctan(e”"’) + 2 arctan(e”)), 


Ripta = Sin(2 arctan(e~) — 2 arctan(e”")), (41) 


fobs 


and if a(t) = e' we have 


R(t) = sin(aresinR +e" —e™), Rpte = Sin(e~™ — e7’), (42) 


3. Asymptotic behaviours of smooth solutions 


In this part, we investigate the asymptotic behaviour of the fields as t > +00. We prove that the 
asymptotic profile U,,(.) := lim,_, ¥(t, .) exists if the scale factor is exponentially increasing. 
In the case of the de Sitter type universe dS,(K), and in the case of the exponentially inflat- 
ing model for which a(t) = e' we are able to compute its expression, and the first terms of the 
asymptotic expansion, in terms of the Cauchy data. Similar results have been obtained in [3] 


for the steady state universe (R, x R$, ds? = dt? — e”dx?) (see also [33]). 


3.1. Asymptotic profile 


We prove the existence of W,, if the scale factor is at least exponentially increasing, i.e. satis- 
fies an assumption of the type: 

a(t) 
a(t)” 


dy>0, 


teR, Vtet* y< (43) 


10 
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Theorem 3.1. We assume (43) is satisfied. Then for any finite energy solution WV of (2), the 
energy (7) is exponentially decreasing 


Vt>t*, E(U,1) <e 2 ME, t”, (44) 
and there exists Y € L(K) such that 
Vrari, VO) —Vllew Se EW, (45) 


Proof. Equations (13) and (43) imply that dE(Ẹ, t)/dt < — 2yE(Ẹ, t), hence (44) is proved. 
To establish the existence of the profile Y, our method relies on the spectral expansion of 
W(t, .) and the energy estimate (or an explicit resolution) for the ODE. More precisely, since 
the Laplace-Beltrami operator Ax on K is a non positive, self-adjoint elliptic operator on a 
compact manifold, its spectrum is a discrete set of eigenvalues —q? < 0. Moreover, by the 
Hilbert-Schmidt theorem there exists an orthonormal basis in L7(K), formed of eigenfunc- 


tions (Uy), C H'(K) associated to q?, i.e. 


—AKY, = Gy. (46) 
Moreover the eigenvalues are explicitly known (see [1, 17, 20, 21]). One has: 


O<q, = 8-1, (47a) 


with 8 € {1, 13,21,25,31, 33, 37, 41, 43, 45, 49, 51, 53, 55,57} U {2n + 1, n> 30}. (47b) 


In particular, since the volume of K is 17/60 (see e.g. [14] p. 5160), we have: 


2V15 
U(X) = ——. (48) 
T 
Any finite energy solution Y(t, X) of equation (2) has an expansion of the form 
W(t, X) = Do ug(tU,(X), teR, XEK, (49) 
q 
where uy is solution of the ordinary differential equation 
/ t) q 
ug 4 at ul Ug = 0. 50 
q a(t) a(t) q ( ) 
For the finite energy solution we have 
2 
2 2 2 2 
EY, 1) => o+ a luo). E O + a? OF < o, 61) 
q q 


and if Y is a smooth solution we have 


dk 
Vk,peEN, oq? ql) <. 
q 
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These series converge locally uniformly with respect to t. Now, multiplying (50) by Us we get 
an energy type equality: 


2 1 2 
roy24 4 2 a(t) Iin) q 2 
|u,(t)| + a(t) |ug(t)| = 2 a(t) [smo a a(t) |ug(t)| } 


We get from hypothese (43) that 


Auot Eok- 2a gc +E uoe 
dr\' 4 a(t) 4 È í N) O S 


hence 


2 2 
u ODP + E lult)? < mal o + a T |ug(t*) e} 


and we deduce that for any t > t* we have: 


q? 2 
Kolce mens ingore J (52) 


Therefore we can introduce 


Ug(Oo) := u,(t*) + ft uj(t)dt, Yao := $ Ug(oo)Y,. 
q 


Now (45) follows from (51) and (52). oO 


We compute W in the next sections. We remark that the profile V,,(.) is more regular 
than the field W(t, .). In short, the theorems 3.3 and 3.5 show that W € H™+! if W(0, .) € H”, 
O;¥(0, .) € H”~ |. This phenomenon has been noticed in the case of the steady-state universe [3]. 


3.2. dS4(K) universe 


We consider the dS,(K) universe for which a(t) = H~'cosh(Ht). We calculate an explicit 
expression of the asymptotic profile and we prove that AxW(t) — za — tanh Ht) '8 V(t) 
tends to zero. 


Theorem 3.2. Given a smooth solution Y of (2), there exist constants A; and A, depending 
only on the initial data such that 


3 
Y(t, X) = |- E Ag + (cosh Hty 243 P, (tanh m| W(X) 
2 


3 3 3 
+ (cosh Ht) 2 Y` lair (tanh Ht) +470? (tanh r| YX), (53) 
q=0 2 2 
2 
q 
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W(t, X) — U(X) = — Z S PATUA — tanh Ht) + o(a = tanh HD ), (55) 
q#0 


1 2o 1 
AxW(t,X) — = — tanh H a(t, X) = aN (1 — tanh Ht% 


+ a Pb gq? + Sajmc0) (1 — tanh Ht) + 0(1 — tanh Ht)2 (56) 
T q#0 


where q and (3 are given by (47a), (47b). The coefficients Aù} are given by 


sin( 2a? +1 ) 
i ie (u0 $ Two) =- r O q=0, (57) 
q+l 


" a 7 
Ag = — |Z uO), A7 = ae) | 0 58 

o=- [SHO p= = a HO. 90 (58) 
where ug(t) := (W(t, .), Uy) Ky 


Proof. The differential equation (50) has the form: 


27492 


ul(t) + 3H tanh(Ht) ul (t) + —t = — 
g(t) (Ht) u,(t) cosh{(H 


Ug(t) = 0. (59) 
For each q we introduce a new function w, defined by 
3 
Uu,(t) := (cosh Ht) 2 w(t). 
Then u, satisfies (59) iff w, is solution of 


Hw} Ct) + (e + >) cosh ~H) — Ang = 0. (60) 


Now we put w(t) = w,(tanh Ht). Thanks to the relation q? + ; =(0 DE + D we deduce 
that w, satisfies (60) iff vw, is solution of the associated Legendre equation 


2 
d— xW) — 2x) + [os 1)— g jw = 0, (61) 
where v = 8 > u= +5 and x € ] — 1, +1[. In the sequel, we use several formulae, all from 


[29]. The first eigenvalue g = 0 is a particular case where v = > Thanks to formula 14.2.7 the 
3 


Wronskian of P20) and P2(x) is not zero. For the other eigenvalues, ( is greater than 13 and 
2 2 


3 3 
is an odd integer. So 07. ı is defined and, thanks to formula 14.2.4, the Wronskian of Pre (x) 
2 


i 
2 


13 
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3 
and Q3 (x) is not zero. So there exist constants A; and A, depending on initial data such that 
2: 


3 aa 3 3 
Wo(x) = Ap P2(x) + Ap P1 2(x), and wx) = AP? ix) + A, Q? (x) for q+ 0. 
2 2 “2 ame 


Then 
3 3 3 
u(t) = (cosh Ht) 2 ai cam Ht) + Ap P, 7(tanh m) , (62) 
2 2 
3 3 3 
u(t) = (cosh Ht) 2 [ar (tanh Ht) + 472? , (tanh n| for q #0. (63) 
-7 -7 


To obtain asymptotic profiles we use hypergeometric representation of the Ferrers func- 
tions (see formulae 14.3.1). We recall that for x€ ] — 1,+1[ 


H 
Po = (1) ro, m1 5-58) 


H 
O (L4x)2 ES w Dr) 1 L) 3 
-(15) 2 Td —p+s) as le 


where (a), is the Pochhammer’s symbol: a(a + 1)(@a+ 2)(a+n— 1)ifn = 1, 2,3,..., and 1 
if n = 0. First we consider the case q # 0, that is v = 8 — 5 = > For t in the neighbourhood of 
+09, the regularity of the previous series at 0 leads to 


3 
1 + tanh Ht i 1 a w+ 1)\(-v) 


3 
2 = 
P;(tanh Ht) = (; ee TCD TO (1 — tanh Ht) 
pet DU HUF D i- tanh HY? FOG — tanh HD») |. 
8r) 
So 
33 3] 1 11 
(cosh Ht) 2P2_ ‚(tanh Ht) = (1 + tanh He} — — (6? — —)—-(1 — tanh Ht) 
= r-i) er) 
(8? — (8-2 
+ (1 — tanh Ht)? + O((1 — tanh H1)3) |. (64) 


TO 
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Likewise for x € ] — 1,+1[ 


3 
3 3 (1+x\i id 1 
Q œ= cos( nf Fw + te) 
6-5 2 sin(=7) 2 —x 272 2 
Twt+ 5) 3 
; l saa n= 
Tp= Gita 22 2 
5 3 
T rw + (i j +00 (w+ =v), 1 at y 
2 Tw- 4) l +x 0 TĒ +s) sth? : 


Then, for f in the neighbourhood of -+00 we have 


-e L z (1 — tanh Ht) + 0((1 — tanh H7)) |, 
ro) 4 205) 


3 = 3 
o? „(tanh Ht) anes ami | 1 
PTa 


2T(8—1)\1+ tanh Ht 


and so 


a 
ITO TO- 


3 5 
(cosh Hty 30? a= = yt tanh At)2 + o(a — tanh Ht). 


(65) 
From (63), (64) and (65) we deduce for t > +00 and q #0 


a 


3 
Uug(t) = -aja + tanh Ht)2 =A; Fe — ha + tanh Ht)⁄2(1 — tanh Ht) 


12T 
+ o(a — tanh H0? ). 


Hence we have got for q = 0 


3 
ug(t) ~ -a f2. and u,(t) + Ay 2 = Pa; Za tanh Ht) + o(a — tanh i), 


(66) 
Using formula 14.10.4 we deduce derivatives of u, from (63). For q = 0 we have 
3 3 
u = = H(cosh Hty2(tanh At)(G — 1) Gow (tanh Ht) + A, p (emh H) 
— H(cosh Ht)-2(8 —1) AtP? (tanh Ht) + AQ? (tanh Ht) 
1 +h 1*1 
3 3 
= H(cosh Ht) 2 (8 — DA; fce HNP? , (tanh Ht) — Par (tanh n| 
-4 1 
3 3 
+ H(cosh Ht) 2 (8 — DA; [amano (tanh Ht) — o? _ a9} (67) 
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And from (67), (64) and (65) we deduce for v = 3 — > a > and t > +00 


u(t) = H(B — paja sna p(t + anh H + (87 D tanh Ht + (G+ B+ b) 
D 


1 fen lp 9 3 S D) 1 
x ro tanh #0) +(e pe- peb- (B+ SVB + VP -p re oon 


+ o(a — tanh #2). 


That is for q #0 


ul(t) = A, 1) (1 + tanh His [o — tanh Hf) — (6? — ba — tanh Ht)? + 0(1 — tanh mo , 


and then 
—2Hq"u,(t) -(d- tanh Ht) 'w/(t) 


ay ae H 3 3 
=q(qit ri WAG + tanh At)2(1 — tanh Ht) + 0(1 — tanh Hr)2. (68) 


It my remains to consider the case q = 0. Actually ug can be a more pD by explaining 


3 
P. We use formula 14.10.4:(1 — efe) (x) = È -> DPIC) + Pia) = $xP2(x). So 
2 


WN! 


P7 has the form C(1 — x?) 3 where C is a constant. Since PO = FÈ, z; 5; 5 = E 
2 2 


we can write Pi) = -2 ad — Py, So with (63) we obtain a new expression of uo: 
2 Gage 
salaj -3AP 2 
uo(t) = —,| — Ag + (cosh Ht) 24y P | “(tanh Ht). 
T 2 
Moreover the hypergeometric representation of Pj? gives us for f at +00 


2 


3 -3 3 
(cosh Ht) 2P, ?(tanh Ht) = (1 — tanh Ht)2 
2: 


1 3 1 
1 Soe ie h At)? + O((1 — tanh H1)3) |, 
TÖ rO Ht) T re: anh Ht) (( anh Ht)°) 
sO 
D o ge | 3 5 
uo(t) = — zæ + Ap T° = tanh Hr)2 + o(a = tanh Ht). (69) 
2 


Then (66) and (69) prove (55). On the other hand we can simplify up since P;? has a simple 
writting (see formula 14.5.18 with 0 € ]0, z[ such that cos 0 = tanh Ht). We have: 
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3 _3 
ug(t) = —3H (cosh Ht) 24 P, (tanh Ht) 
2 


3 3 
= —3H (cosh Ht)? ——Aa (sin 02 =H [2 (cosh Ht) 3Ap 
TC) T 


3 3 
= -H 2 ati + tanh Ht)2(1 — tanh Ht X2, 
T 


So 


3 1 
Ao (1 + tanh Ht)2(1 — tanh Hr)2. 


1 1 
—— (1 — tanh Ht) luh (t) = 
2H 0 Jm 


Finally to compute AZ we remark that since 8 = Jq? +1 is an odd integer number, we 


have 
23 1 z3) 2 
P1200) == /=, |P17|(@ =-,|=, (70) 
z 2\2 A T 
3 7 rT 3 y 
P? ,(0)=—,/—Gsin| B—}], | P2 0) = 0, 
210 2 osin(93) | a ) (71) 
3 a | T 3 T 
2 (0) = 0, |Q? 0) = —,/—q° sin| G— |. 
EQ le ) iF 4 sin( 3) (72) 
The expressions of the coefficients follows from (53) by elementary computations. The proof 
is complete. O 


We deduce from this theorem that any smooth solution Y can be expressed as 
Vt, X) = VH, X) + V, X), 
where Ų+ and Y~ are smooth solutions satisfying 


wr, X) = U(X) + Oe, (73) 


wt, X) = o(e2""). (74) 


It is sufficient to define 


3 
wer, X) := — E Aj W(X) + (cosh Ht)"3 5y AP? (tanh HAWX), (75) 
q0 2 


-3 3 
Y(t, X) := (cosh Ht? jer, *(tanh HN Wo(X) + > A, Q?_ (tanh HY |. (76) 
2; qz#0 72 


We express the initial data for Y= in terms of those of Y, and we estimate the norms of Ya. 
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Theorem 3.3. The smooth solutions Y, U~ and the asymptotic profile V,, satisfy : 


U+(0,X) = Y0, X) + T (awo, ), Yo)zaot(X), 3V0, X) = 0, (77) 
U-(0,X) = Faw, ), Vo) 24) %AX), 3V0, X) = 3V0, X), (78) 
T= re ar sin = De FT jro, ), (79) 
Zll = ITO, Vlla» (80) 


[| Vic Wocllamacy = | VkY +0, Jlla-a9 = IVKUO, Jla- m= 0,1. (81) 


In (79), the pseudodifferential operator f (Ax) is defined by the usual functional calculus 
of the spectral theory: f(Ax)® is defined by (f(Ak)®, Yq) pK) = f(—q’)(, Wy) 2K). We note 


that (47a) assures that sin( 2a? +1 ) €{+1, -1}. 
Proof. We obtain (77) and (78) from (70)-(72). We get (79) by (77) and the equalities 


w+0,x) = —,|2 Essiz Jaz, = ed sin( = JD + TNs. 
T 
q 


Finally since sin(23) = +1, (80) and (81) follow from (77) and (79). ðO 


2 


We remark that Y- = 0 iff &¥(0, .) = 0, and the asymptotic profile Y depends only on 
W*(0,.), therefore it is determined by (0, .) and h aV (0, .). Moreover the map 


S: VHO, .)— T 


is an isomorphism from H- !(K) onto L(K) and from Z?(K) onto H'(K). Since S is a pseudo- 
differential operator given by (— Akp + 172 TEN —Ak +1 ) the computation of the 


asymptotic profile using this formula is a very hard task. In contrast, the numerical method for 
the computation of the time dependent field Y(t, .) presented in the next part, allows one to find 
W(X) = lim,., +o V(t, X) with a good accuracy. 


3.3. Inflating universe 


We consider the toy model of an exponentially inflating universe with a(t) = e’. We calculate 
the asymptotic profile of a wave W and we show that Ax W(t) — e”3,Ņ(t) tends to zero. 


Theorem 3.4. Given a smooth solution Y of (2), there exist constants A; and A, depending 
only on initial data such that 


W(t, X) = [Ape + At] WX) + e7? > [a 3(ge) + A7 Y3 (aem) WX), (82) 
q=0 2 2 
U(X) = Ap U(X) — = > q 2ASY,(X), (83) 
q0 
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Wt, X) — U(X) = CRE 25 ary, | + 0(e~*), (84) 
2\a0 q#0 

A(t, X) — e”a Wt, X) = TES Jald Z alaro] + 0(e~7/) (85) 
qs 


where the set of q is given by (47a), (47b). The coefficients Az are given by: 


Aj = FGO + 0), Ap = — > G sing )uj(0) + 4 $004 Ho} (86) 
ga 


Ay = O; A = as G £08) (0) + sca ak LF io} (87) 
4 2 E 


where u,(t) := (W(t, .), Vy) pay 


Proof. If a(t) = e’, equation (50) is just 
ug (t) + 3uj(t) + gre "ug = 0. (88) 
For each q we introduce a new function w, defined by 
Uu,{t) = e 2 W(t). 
Then uy satisfies (88) iff w, is solution of 


wi (t) +(e 26 —2t _ = hi —0. (89) 


q = 0 is a particular case. wo can be written wo(t) = Age? + Agetz where Aj and Ap are 
constants depending on initial data. So 


u(t) = Ape" + Ag. (90) 


Then we consider q = 0 and we put w(t) = w,(qe~"). wg satisfies (89) iff vw is solution of the 
Bessel equation 


2 
x 2 (x) + xwi(x) + [e — (ż) jaw =0. (91) 


For x belonging to ]0, +-oo[, W can be written W(x) = Ay Y3(x) + A, J3(x) where A, and Ay 
2 2 


are constants depending on initial data (see 10.2 of [29]). So we have 


ug(t) =e ue gJ3(ge) +e Ay Yaqo- 5, 
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According to 10.6.2 ([29]) the derivatives of J3 and Y3 satisfy C(x) = Ci(x) — +C3(x). We 
obtain 7 2 2 7 7 


3,) 3 aa 
u(t) = —e JEZ N+ qe Aad ae + 54a rage ") + ge'A; re ) 


= —34 = —t mer + —t 
=-e 2qA,Ji(ge y-e 2 q A; Yi(ge ). 
2 2 


In fact, by the use of some formulae in 10.47ii and 10.491, these Bessel functions have simple 
expressions: 


Ji(x) = A Jo(x) = sinx, Y\(x)= 2x Yo(x) = — 7 cos x, 
2 \ r \ mx 2 \ r \ 1x 
2% 2x | sinx cosx 2x 2x cos sin x 
x(x) = |= A@ =| | : ] 3 = | ya) =|] | * | 
2 T T x x 2 T T X x 


and we deduce (86), (87) by elementary calculations. To investigate the behaviour of ug at +00 


we use the following asymptotics at 0 of Y3, J3, Yı and J1. 
2 2 2 


2 ļ| 4 15 1 9 13 
Ji (x) = E x2 + 24a), 
2 T 6 120 
3 15 
Yı(x) = : | x 3+ La Li l vt + 00], 
2 T 2 24 720 
| 3 
J3(x) = 2 Ee l x + : x2 +062)], 
2 m L3 30 840 
2 3 Lr. 15 1 9 13 
Y3(x) = | X 2 X2 + —x2 x2 + ow] ; 
2 T 2 8 144 
All these relations show that for any q = 0 we have at +00 


u(t) = E e73'A7 [5e] $ E erat |-aey2 - Lq} + rae? + oe*h 
T 3 T 2 8 


2 ei l 1 1 3 
EI q 2A; Lgzaje + daare] + oe, (92) 


and 


5 2 1 5 2 yy 1 
u(t) = —e Hal? 4, (qe +e val Ai (qe) 2 + 0(e-*) (93) 


3 
= 2 [24e 1 — q747e * 1+ o=. (94) 
T 


(84) is proved thanks to (90) and (92). Moreover for all q # 0 (94) gives 
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i 3 
-gut — etua = ,|2 grat — 2 224; — qiAze [+o 21 
T Kig 
DOEDE m 
= — qe ‘+ Q(e~~), 


AKW(t,X) — e OME, X) = eu NWX) + SO (qult) — e”u (MUX) 
qz0 


Then, thanks to (90), we deduce 


= (3Ay Uo(X)e“ | => auo) +0. 
qz#0 


This achieves the proof. O 


This theorem shows that the situation for the inflating universe is similar to those of the de 
Sitter space. We can write again 


Wt, X) = V(t, X) + V(t, X), 


where Y+ and Y~ are smooth solutions satisfying 


VH, X) = U(X) + Oe”), (95) 
w(t, X) = O(e*. (96) 
It is sufficient to introduce 
3 
VH, X) := Ap Wo(X) + e77 X A} Y3 (qe YX), (97) 
q#0 2 
3 
W(t, X) := Ape WX) + e77 X A7 J3 (qe DNA). (98) 
q#0 2 


The coefficients Az are calculated from the initial data of Y at t = 0, but we can also deduce 


them from (Y+(0, .), ôV +(0, .)). Y*(0, .) characterizes Aj and the family Ay Y3(q). O,Y*(0, .) 
2 
characterizes the sequence ahina + arso} Since |¥3(q)| + | Y4(q)| = 0, the whole 
2 7 2 3 
family Aj is determined. To estimate the norms of the asymptotic profile Yx, we introduce 
generalized energies at time zero of the solutions Y of (2) by putting 


Em (Y, 0) := ||VK Y0, Vlla + WOO, Vra, m € Z. (99) 


Theorem 3.5. There exists C, C' > 0 such that for any smooth solutions we have 


X+, )lfza + CW E208, 0) < IY) < WTO, Iza + CE-2(*, 0), 
(100) 


CE_\(U*,0) < || Vrtal) < CE1(U*, 0), (101) 
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CE (U*, 0) < JAK Yll?) < CEo(Y*, 0) < C'E, 0). (102) 
Proof. We expand the norms in terms of series: 


2 2 

2 = 2 = 

Wolleay= [Ad P+— > q JAF g IV Yll) = = q AZP, 
T q#+0 T q#0 


2 
AkVcllzay= = >> qA} P, 
T q#0 


2 _ . 
W*(0, llarg = Ag P + Xa + D"q*|cos q — gsing?|Az/, 
qz0 


2 ; cos 
Ent, 0 = = Sq? + vfa — sin(2q) + a); P, 
q#0 


co: 


and we use the inequalities c~'g<q—sin(2q)+ 7 Sed for some c>0, and 


|A} P S qalu (OP + q'u OP. oO 


4. Numerical resolution 


We develop in this part an accurate scheme of computation of the wave propagating on K in 
both cases a(t) = H~! cosh(Hft) and a(t) = e'. The computations are performed on F,. The pre- 
cise computation during a long time is a hard task because of the exponential behaviour of the 
scale factor. To overcome this difficulty we use finite elements of second order. We carefully 
check that the numerical diffusion is very weak by computing waves with a future horizon, 
and we compare the numerical localization of this horizon with its theoretical value given by 
(36). The validity of our approach is established by comparison with explicit solutions, and by 
numerically testing the asymptotic results of the previous section. In particular the asymptot- 
ics (56) and (85) show that 


A(t, .) — ADO, .) = Oe), t> +00, (103) 


for some a > 0 where A is an exponentially increasing function. Since 0,W(t, .) tends exponen- 
tially to zero, the numerical checking of this property is challenging: it is an excellent test of 
the robusteness and the accuracy of the scheme. 

We solve the variational problem (30) using the usual way (see e.g. [10]). We take a family 
Vi, O < h < ho, of finite dimensional vector subspaces of W'(F,). We assume that 


Tv, — wl 
zis Va = Ww (F). 
We choose sequences uo,h, U1, € Vp such that 
uon > Uo in WF), wpn > u in D(F). 


We consider the solution u, € C®(R;; Vp) of 
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2 
YEW 05$ f. O- IXPY zul, 4004X 


add 2l 
“a a J, OIP tut. X06 (20K 


J, (1—|XP)-2Vup(t,X) - Vot, X) dX 


Fo 
pees _lyl2\-s(x. 
25 f.o |X|) 2(X Vult, X))(X - Vo, (t, X)) dX, 


satisfying up,(t*, .) = uon(.), Orun(t*, .) = u,a(.). We shall put an initial data at t= t* with £* 
large enough, rather than at t = 0, in order to achieve faster limit states. We consider a basis 


h : so 
(e a <i<m of Vp and we expand up on this basis: 


Nn 

h h 

up(t) = 5 uj(t)e;. 
jal 


We introduce 


.— t¢,,h ,,h h _ = _ 
U(t) amn: (uj, U, tt, un,)» M —= (Mi) <ij<my ) = Dip) <ij<my K = (Kii <ijcn, 
where 


B 1 Rescate 
ij =Í ep O 


1 h h h h h h 
y= | —— 0O) + 8l lK) + O,e"(X)Oye"(X)) aX, 
i I, J1-|XP , ai , 


1 


(xO,e4(X) + yOye(X) + 20-e)(X)) dX. 


(xO,e}(X) + yOye(X) + 20-e7(X)) 


Then the variational formulation is equivalent to 


a Oy 1 

M| U” + 3—— + K+D)U=0. 104 
l a(t) v) Po T i 

Like in [5] we can prove that the solution u, obtained by solving (104) tends to the solution u 

of (30) as h > 0. This differential system is solved very simply by iteration by solving 


a'(t* + nAT) r v| yrti_ | 
2 


0= M U”"+l— 2U" + un} +3 

( ) a(t* + nAT) 

1 

+ (ATY————_——_{K + D)u" 
of) TP + nAT) ) 
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with an initial data at t*. The Lax—Richtmyer theorem assures that the approximate solu- 
tion given by this scheme tends to up when the time step AT > 0. Moreover we compute 
Eq(t* + nAT), an approximation at time t* + nAT of the energy E(t* + nAT);: 


(K | D) ur! U”). 


(105) 
We use the same mesh as in our previous article [5]. But we construct the finite element 
spaces V, of P, type instead of P, type. We perform this choice to gain a better accuracy 
that is necessary to overcome the difficulty linked to the exponential behaviour of a(t). We 
take into account the boundary condition (24) in the definition of the finite elements, so that 
V, C W(F,). We note T, the set of all tetrahedra of the mesh, Fn the set Urez, T, and P(T) 
the set of second degree polynomial functions on T. We choose Pz Lagrange finite elements, 
so our nodes are among vertices and middle of edges of T€ 7}. In our construction of the 
mesh we took care that vertices belonging to 0(F, n) belong to O(F,) and that for any vertex 
S of TE Thn, S’ is also a vertex of a tetrahedron of T} if S’~ S (see [5]). It is therefore easy to 
take into account the boundary condition for vertices of T € 7}. Moreover middle of edges of 
peripheral tetrahedra are not on the border of F, thanks to the convexity of edges and faces of 
Fy, so they are only equivalent to themselves. Then we introduce: 


n_ ypyn-l n_ ypyn-l 
ae +nary= (M R ) l 


At , At a(t* +nAT) 


Vr := {v : Fun R, v E CAF, a), VTE Th, vir € PAT), M ~ M' => vM) = v(M’)}. 


dim V, is equal to the number of equivalence classes of vertices plus middle of edges of tet- 
rahedra of Tp. If j is the number of a node associated to a vertex M; or to the middle A; of an 


edge, we construct a basis (e), of V, by: 


<I<Nn 


(i) If j is associated to a node M; that does not belong to OF, : e(M;) = by, e(Ai) = 0. 
There are n,; functions of this kind. 
(ii) If j is associated to a node M; that is a vertex of Fy: 
h 1 ifM ~M, p 
e(M;) = ,e(A;) = 0. 
ia f otherwise iA) 
There are five functions of this kind. 
(iii) If j is associated to a node M; that belongs to a face of OF, and not to an edge: 
ei(M;) en — Py 
0 otherwise 
There are 6 x n, functions of this kind. 
(iv) If j is associated to a node M; that belongs to an edge of a face of OF, and is not a vertex 
1 ifM;=P, 
0 otherwise 
There are 10 x n,e functions of this last kind. 
(v) If j is associated to a node A; that is a middle of an edge: ei(Ai) = ĝi ei(M;) = 0. 
There are as many functions of this kind as different edges in the mesh. 


(Ai) = 0. 


of F;: ei(M;) = { , eN(Ai) = 0. 


with nye the number of mesh vertices on an edge of a face that are not a vertex of F,, my the 
number of mesh vertices on a face that are not on an edge and n,; the number of mesh vertices 


O 
in F,. For a given mesh, the number of nodes is important. For example if we have a mesh with 
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Time 


Figure 3. Solution for Init), on the left at t* = 3.5, on the right at fo, = 28.7. 


89 vertices on each edge of OF,, we have 730309 nodes for P, type finite elements method, 
and 6377052 nodes for P, type method. Matrices are created with a 31 points tetrahedral 
quadrature formula which is of degree 7 ([{18]). 

With our choice of finite elements, the Laplacian is constant in a tetrahedron 7; in the fol- 
lowing it will be noted Azup\r(t, Gr) where Gr is the gravity center of tetrahedron T. To test 
(103), given by (56) and (85), we have chosen to calculate Norm defined respectively by: 

1 
2 


2 
Norm := (Arn Gr) - sl — tanh Ht) !O,up(t, cn) dX 


pa aj 1—|X/? 


and 


Nie 


(Agupyr(t, Gr) — eO,up(t, Gr) aX 


1 
Norm := Í —— 
p2 T y 1- |x}? 
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3 


Figure 4. Solution for Init, on the left at t* = 1.5, on the right at fons = 19.7. 
4.1. Numerical results for the exponentially inflating universe 


For this model the scale factor is a(t) = e'. First we check the efficiency of our code by the 
computation of the trivial solution Y(t, X) = 2e~*”. The figure 2 shows that the result obtained 
by using our code matches with the exact solution e~*” depicted with MAPLE. The results 
agree at +1075. 

We now test the accuracy of our code by the numerical localization of a future horizon and 
its robustness by checking the asymptotic behaviour (103). We choose different initial data 
and a small time step: 0.000 15. We present our results for two initial data at time t* denoted 
Init, and Init, of the form: 


#(X,X0) 
ug(X) = Ape#%.Xo)—R> for d(X, Xo) < Ro, and uo(X) = 0 for d(X, Xo) > Ro. (106) 


We take Ap = 100 and Xo = 0. For Init; we have Rọ = 0.1 and for Initz, Rọ = 0.05. In order 
to simplify we choose O,u(t*, .) = u,(.) = 0. Figures 3 and 4 present the initial data and the 
limit state on the slice z = 0 of F, observed at time tops > f* that we choose large enough in 
order to that the asymptotic state is reached. Figure 5 shows the time evolution of the solution 
associated to Init; for some points P of F, pointed out in figure 3, and figure 6 shows the time 
evolution of the Init, for some point PP of F, pointed out in figure 4. In the following, MO 
denotes a node of the mesh that is the closest to the center 0 of F,,. 

There is a future horizon and R, = sin(arcsin 0.1 + e~3>) œ 0.13. 

In this case we do have R, = sin(arcsin 0.05 + e~!°) ~ 0.27. 

It clearly appears that a limit state is reached very quickly, as soon as t = 8 in both cases. 
We also note the existence of a future horizon. In tables 1 and 2 we can read the distance 
between previous points M(xy, yy. zm) of F, and its center. The agreement between the value 
of R, given by (42), (41) is quite good, but it is difficult to have a precise numerical estimation 
of this horizon; for example P2 and P3 belong to a same tetrahedron, so the theoretical horizon 
passes throught this tetrahedron. For Init, and t* = 3.5 it is found that outside the ball of radius 
0.129 the solution is always less than 0.01, and outside the ball of radius 0.142 the solution is 
always less than 0.001. Knowing that the limit state reaches 90 we conclude that the numerical 
estimation of R, is satisfactory. 
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Figure 5. The solution u(t, X) at MO, P1, P2, P3 with Init, and ż* = 3.5. 
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Figure 6. The solution u(t, X) at MO, PP1, PP2, PP3 with Init, and t* = 1.5. 
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Table 1. Init; and t* = 3.5, R, ~ 0.1299. 


Point dys. M) u(15, M) 
MO 0.000 08 90.45 
Pl 0.064 < R 37.47 
P2 0.126 < Rn 0.008 


P3 0.131 >R, 0.0016 


Table 2. Init and r* = 1.5, R, ~ 0.2697. 


Point d(0_3,M) u(15, M) 
MO 0.00008 < R 0.0032 
PPI R < 0.218 < R, 0.51 
PP2 R < 0.245 < R, 0.22 


PP3 R < 0.270 S Rn —0.0017 


For large time, the time derivative becomes smaller and smaller, so the energy continues 
decreasing and there is no longer evolution even to t = 30. Ey defined by (105) is decreasing. 
For example for Init, and t* = 1.5 we have obtained: 


Timet E4) Timet E4) Timet EÐ Timet £,(t) Timet E(t) 
t* = 1.5 217.16 3 0.56 48 4.493 11 2510-9 19.5 107! 


The test of the asymptotic behaviour (85) is presented in figure 7 where we note that Norm is 
rapidly decreasing until t = 10. We emphasize that the exponential growth of the coefficients 
makes the computations very difficult, and this test shows the robustness of our numerical 
method. 


4.2. Numerical results for dS4(K) 


We assume the Hubble constant is equal to 1 hence the scale factor is a(t) = cosh(ft). First 
we test the accuracy of our scheme in the trivial case of the constant solution W(t, X) = 2. 
Taking (0, X) = 2, W(At, X) = 2 with Ar = 1.5 1074, our scheme gives a constant solution 
Wt, X) = 2+ 107! for t € [0, 3.5] and the energy is around 107". 

On the other hand we test our scheme by comparing it with two other numerical methods 
for computing the solution 


3 
W(t, X) = Ao (cosh NZP] 2(tanh 4). (107) 
2 


In figure 8 we present the results computed by: 


(i) our finite elements scheme with initial data 


(0, X) = ` — , W(At, X) = WO, X) — {2 At, At=1510-4 (108) 
win 


(ii) a Runge—Kunta scheme for the ordinary differential equation (59) with the same data 
(108). 
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Figure 7. Validation of asymptotic behaviour (103): Norm with respect to time. 
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Figure 8. The black curve is the superposition of the three colored curves. 


Gii) a symbolic calculation of (cosh 1) 3P;>(tanh t) by using MAPLE. 

2 
The three curves perfectly matche and we can see only a unique graph. In fact the numerical 
values are the same at +5 10~°. 

In the sequel we take initial data Init and O,u(t*, .) = 0, for various f*. This is the initial 
data presented in the first picture of figure 4. First we check the accuracy of the numerical 
localization of a future horizon. We choose t* = 3.5. Then R} ~ 0.11. Results are presented 
in figures 9 and 10. Figure 9 shows the initial data at time f* = 3.5 and at time tobs = 20.08 
for which the asympotic state is reached. The pictures are on the slice z = 0 of the dodecahe- 
dron. There is a future horizon and the support of the field is strictly included in K. Figure 10 
presents the time evolution of the solution at some points MM of F, depicted in figure 9. In 
tables 3 and 4 we can read the distance between previous points M(xy, yy, zm) of F, and its 
center, and all these numerical results are in good agreement with the theoretical results on the 
future horizon (36) and the asymptotic profile. 

There is a horizon and R, = sin(arcsin 0.05 — 2 arctan(e*>) + m) œ 0.11. 

As for a(t) = e' we see that the limit state is quickly reached: as soon as t = 9 in both cases. 
E, defined by (105) is decreasing as we can see for example for Init, and t* = 1.5: 


Timet £,(t) Timet E4 Timet Eð Timet Ea) Timet Eð 
t* = 1.5 4361 3.4 0.49 5 1.21073 11.5 10° 19.6 9 107" 
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Figure 9. Solution for Inito(Ro = 0.05), on the left at t* = 3.5, on the right at 
fobs = 20.08. 
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Figure 10. Solution u(t, X) with initial data Initz at t* = 3.5. 


Finally we test the robustness of our code by the hard computation of (103). We consider 
the solution with initial data given at f* = 1.5 and f* = 3.5. In all cases the quantity Norm that 
tests (56) quickly tends to zero. On figure 11 we remark that for t* = 1.5 Norm is decreasing 
except when the support of the solution reaches the border of F, at t ~ 2.3 (we detail the pecu- 
liar properties of this solution below). We conclude that as for a(t) = e’, the tests of (56) by the 
computation of Norm prove that our numerical scheme is accurate and robust. 

We end this part by the most interesting case: we present three solutions with ‘the circles- 
in-the-sky’. Firstly we consider the previous solution with data W(¢*, .) described by Initz and 
O,W(t*, .) = Ogivenatt* = 1.5. Then R; is equal tosin(arcsin 0.05 — 2 arctan(e!>) + m) ~ 0.469 
which is greater than sin(dmax) = 0.378. Therefore the support of the solution is the whole F, 
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Table 3. Init and t* = 1.5, no horizon. 


Point d(Ogs. M) u(15, M) 
Mo 0.000 08 —0.0005 
MI 0.270 0.415 
M2 0.182 0.004 


M3 0.354 0.013 


Table 4. Init? and t* = 3.5, R, = 0.11. 


Point d(Oy3, M) u(15, M) 
MO 0.00008 —0.016 
MM1 0.0544 7.41 
MM2 0.1004 0.131 


MM3 0.1055 —0.008 


for t large enough and the constraint (40) is fulfilled: we shall be able to see multiple copies 
of the asymptotic profile. 

Figure 12 shows the solution at time fobs = 15.06 on the boundary of the dodecahedron, 
OF,, and in its interior, on the slice z = 0. In fact, at this time the asymptotic state is reached. 
We can see this fast convergence on figure 13 that shows the time evolution of the solution 
at some points M of F, depicted in figure 12. Since the whole universe is included inside the 
horizon sphere, it is interesting to present the pull-back of the solution in the universal cover- 
ing S of K. Figure 14 shows the solution on a partial tiling of the unit ball B’ (see figure 1 and 
the explanations after formula equation (14). 

We now show the application of our numerical method to the determination of the temper- 
ature fluctuation in the deep sky. We consider an observer located at time tọ at the centre of the 
Poincaré dodecahedron, x = y = z = 0, that corresponds to the north pole of S°, defined by 
(x9, x!, x2, x3) = (1,0, 0,0), or in spherical coordinates, y = 0. We assume the initial pertuba- 
tion of the metric is given by Init, at t* = 1.5 and tọ is large enough for the asymptotic state to 
be reached. Then according to the Sachs—Wolfe formula (4), the temperature fluctuation seen 
by this observer in the direction (0, p) € S? is given by 


6T(6,~) _ 1 
= —aa0 a —W,, k 0, k 
- 3 (x, 9, p) 


where the angle x is determined by X = no — Mı, that is the conformal time between the time 
of recombination tı, and to (t* < tj; < to). Its value is given by equation (41), 


x = 2 arctan(e”) — 2 arctan(e“), 


and the sky seen by the observer corresponds to the horizon sphere centered at the origin, 
included in B’, with the radius R,,.,, = sin x. Figure 15 shows the deep sky for Ry, = 0.5 and 
Rit) = 0.6 (we recall that by (31) the radius of the smallest ball containing K is Rmax % 0.378). 
We remark that the initial perturbation of the metric leads to a temperature fluctuation that 
shows spherical correlations along six pairs of antipodal matched circles (the famous circles- 
in-the-sky) that are a signature of a complex topology with a positive curvature. In contrast, a 
complex topology and a negative curvature lead to chaotic temperature fluctuations [4]. 

We now present a more complex solution associated to an initial data u(t*,.) = Init, 
O,u(t*, .) = 0 given at time f“ = 2. Initz is composed of Init, to which we have added two another 
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Figure 12. Solution observed at fobs = 15.06 for Init) and t* = 1.5. Left: on the 
boundary of the dodecahedron. Right: on the slice at z = 0. 


fields of the same type (see (106)): one with Xo = (0.1, 0.1, 0.1) and Ro = 0.1, the other with 
Xo = (—0.15, 0, —0.1) and Ro = 0.05. We take Ap = 100 again. Figure 16 shows Init; and the 
asymptotic profile on the two-plane 10x + Sy — 15z = 0. The left figure 17 shows the pull-back 
of Y on a part of the universal covering S? of K, included in B’ (see (15)), composed by F,, 
P(g(F)) with i = 1,..., 12, p(g, © 96(F)) with i = 1,...,5, p(g; © g,(F)) with i = 1, 2, 4, 5, 6 (see 
(18) and the appendix). On the right the deep sky for x ~ 0.582 is depicted. Similarly figure 18 
presents the case of a large conformal time. On the left, the part of the universal covering 


of K is included in the second unit ball B” associated to x € [5> r]. This set is composed by 
P(-IF), p(g(—1F)) with i = 1, ..., 12, p(g, ° g6(—1F)) with i = 1, ..., 5, p(&; ° g3(—1F)) with 
i= 1, 2, 4,5, 6. On the right the deep sky for x ~ 2,498. We note the correlations along six 
pairs of circles again. 

We achieve this part with the more realistic case of a random initial fluctuation Init, com- 
posed with 100 localized fluctuations of type (106). The radius Ro, the centers Xo, and the 
amplitudes Ap are randomly choosen such that — 100 < Ag < 100, Ro < 0.1. The left figure 19 
shows Init, on a cut of the dodecahedron along three planes. We compute the solution for 
u(t*, .) = Init4, O,u(t*, .) = 0 given at time f* = 2. The right figure 19 presents two antipodal 
views of the deep sky for x = 7/2. We can guess the correlations of the fluctuations along six 
pairs of matched antipodal circles. 
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Figure 13. Solution u(t, X) with initial data Init at t* = 1.5. 
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Figure 14. Pull-back of the solution at fobs = 15.06 in a part of the universal covering 
S? of K. 
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Figure 15. Deep skies: on the left x = arcsin(0.5)~0.523, on the right 
x = arcsin(0.6) = 0.643. 


Figure 16. Solution for Init; on the cut 10x + 5y — 15z = 0: on the left at f* = 2, on the 
right, the asymptotic profile. 


Figure 17. On the left the asymptotic state in a part of the universal covering S* of K 
around (1, 0, 0, 0). On the right the deep sky for x = arcsin(0.55) ~ 0.582. 
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Figure 18. On the left the asymptotic state in a part of the universal covering S? of K. 
On the right the deep sky for x = m — arcsin(0.6) ~ 2.498. 
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Figure 19. An initial random fluctuation and two antipodal views of the deep sky for 
xX = 1/2. 


5. Conclusion 


In this paper we have considered models of dynamical spacetimes for which the spatial sec- 
tion is the Poincaré dodecahedron. We have solved the Cauchy problem for the D’ Alembertian 
and expressed the initial value problem in terms of variational formulation put on the funda- 
mental domain. We have established the existence of an asymptotic state Y, at the time infin- 
ity for the finite energy waves: 

Py := lim Y(t). 

f=+00 

If W is the perturbation of the metric, this state describes the temperature fluctuations in the 
deep sky thanks to the (ordinary) Sachs—Wolfe formula 

6T 1 


— x -i 


T 3 
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We have obtained the analytic expression of this limit state for two important cases of acceler- 
ating universes associated to the scale factors H~! cosh(Ht) and e’. It turns out that Y is given 
by a pseudodifferential operator acting on the inital data of the perturbation, ((0), 0,Y(0)) that 
is not convenient for the computation. Therefore we have performed a numerical method for 
solving the time evolution problem. We also emphasize that this step is necessary if we want 
to use the integrated Sachs—Wolfe formula involving W(t) for any time ¢ since the last time of 
recombination. We employ the P, type finite elements for the sake of accuracy. The numerical 
experiments show a good agreement with the theoretical results on the asymptotic behaviours. 
That proves the robustness and the accuracy of our method as well as its efficiency to compute 
the deep sky. Its drawback is the lenght of the computation, usually we need one week to go 
from ¢ to t+ 1.2 and this time of calculation would explode for a very refined mesh. We can 
hope to be able to use a parallel computing together with a domain decomposition on a sig- 
nificantly refined mesh. Then we shall be able to treat the heavy computations associated to 
random initial data and the integrated Sachs—Wolfe formula. 
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Appendix 


This appendix is devoted to a brief description of F and F,. For more details on the construc- 
tion of these domains we refer to [5]. 

We warn that figure Al is not a faithful representation of F, as, for convenience, we have 
replaced the true curved faces of F, by the plane faces F j that are the projection on R? of the 
sets of all barycenters in R4 of the vertices S£, ar s of the faces F; of F. 

We need the four basic quaternions {1,i, j,k} of the set H of all quaternions to define the 
twelve Cliffort translations g; involved in the construction of F: 


1 11. 1 1 1. 11 1 ci. 1 

= —ol + i+ —j, := —ol+ k, := —ol+ j k, 
8 5 te 2? 273 2 20 eae) 20) 2 
1 11. 1 1 11. 1 1 1. 11 

P pee es | ee gh i a fe 
an: 20 2P 8797 2a 3 Bo a” Be 


They are such that: Vie {1,...,6}, g;(FE) = Fi+6. The inverses of these six first translations 
are the six other translations: 


1 11 1 1 1 11 

:= (¢,) ! = —o1 — —-i- >j, = (y! = —ol — —i+ ——k, 
= (B= 5 23o gh 8 = 5 2 20 
1 11 1 1 11 1 

=f ° r —1 i i 

= ==¢1 = jk = = 701+ ——-i- —j, 

go = (83) 5 ne ee (84) J ao) oO 
1 11 1 1 1 11 

= ~! = øl j k, = t= øl i . 
811 = (85) 7 7 z! 7 812 = (86) 7 2 20 
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Figure A1. Diagram of vertices S;, edges and faces of F,. 


(i) Coordinates of the vertices S; of F: 


1 1 1 1 1 1 
S a, ’ > } S = sale 1, aR o} 
' ll oo o ? 2/2, o 
1 1 1 1 1 1 1 1 
S o°, ; ; } S4 = (2. ; i ) 
‘ zl o OO i 2/2 o o o 
1 1 1 1 1 1 
Ss = o”,0,—1, } S6 = (<°. ; ) 
? -zl o? j SyF ooo 
1 1 1 1 
S7 = —— o?,— 5.0.1}, Sg = zhao, L5) 
7 2 ( o? 5 2/2 o? 
So = : (2 : > : ’ : } Sio = sal” Zo i} 
2/2 ooo 242 o 
1 1 1 1 
S = — a, 0, 1, -4} S = sal —-1, $ o} 
5 zal Pr H a 2 
1 1 1 1 1 
S3 = o, 0, i} Sig = (2. ; } 
i mal 2 a 2/2 o oo 
1 1 1 1 1 1 
Sis = a”, —,0, i) Sig = (0°. : , } 
: =a o? ee oo o 
1 1 1 1 1 1 
Si7= eS ) Sig= (0°. i o); 
i mal oo 8 I2 a 
1 1 1 1 
AY = — 2 : ——,0 S = 2, 0, l, 
i zal’ o a TA > 


A Bachelot-Motet and A Bachelot 


(ii) Images by the Clifford translation g; of the face F; of F and of its edges. 


gı maps F; to F}, and we have for the vertices and the edges: 


&(S3) = Se, (Sis) = Sg, 9(Si6) = Sin 085) = Si 8(S20) = Sd, 
8 ($3818) = SoSg, 80518516) = SgSi1, 3CS1655) = SipSi7, g0S5S20) = S1782, 


82053) = S286. 
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g2 maps F to Fg, and we have for the vertices and the edges: 

&(Sis) = Sis,  8(Si2) = Sir, 8S9) = S2, 8S7) = Sio, 8(S3) = S4, 
go SisSi2) = Sis8i7, 81752) = SeSi» Bo(S0S7) = S2819, 85783) = SoS, 
9551s) = SySis. 

g3 maps F; to Fo, and we have for the vertices and the edges: 

83(S3) = Si, 83(S7) = Sin 83(Sio) = Siz, -83(Si4) = Sis, 83(S20) = Sis, 
BEN) = SiS, BES) = SaS, g4(SioSia) = SirSis, 

g3(Si4S20) = S15S13; 83(S2083) = Sias). 

g4 maps F; to Fo, and we have for the vertices and the edges: 

84(S20) = So,  84(Si4) = Sg, 4(Si9) = Sir, 84(S4) = Si, 84(Ss) = Sia, 
B4(S20814) = SoSs, B4(SiaSi9) = SeSi 84(Sio84) = SuuSi, g4(SuS3) = SiS» 
84( 55520) = $1289. 

gs maps Fs to F11, and we have for the vertices and the edges: 

85(S5) = Sio,  85(S4) = S6, 85(Sis) = Sg, 85(S13) = So, 85(Si6) = Sz, 
85( 5552) = SioS6, 5(S45i5) = SoSs, 95(SisS13) = SeS» gs(SiaSi6) = SoS, 
95(SioSs) = SSi: 

gs maps F¢ to F2, and we have for the vertices and the edges: 

gelS16) = Sio,  6lS13) = S2,  86(S1) = Sos B6(S12) = Sios Be (Sig) = Sia, 
Bo(Si6513) = 1S2, 86135) = S286, g6(SiSi2) = SeSi0, 


g6(Si2Sis) = SioSia, g6(SisSie) = SiaSio. 


Then we define the equivalence classes on F by: 


qEF 34={4), 

AGI‘), QERNKOK => 4 = {q, 8:4), g4), &(Q)}, 

(AG), qERNK) and Yki, j, qg F) >4= (4,8(9),8(Q)}, 
Gi, qEF) and Wii, q¢F) >4={4,8{q)}. 


Here the integers i, j, k belong to {1, ..., 12}, and F; denotes a face of F. 
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